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1. Introduction 



Quantum Field Theories with boundary conditions Q |§ § || § have interesting 
applications to quite a large class of physical phenomena, among them the Kondo problem 
, absorption of polymers on a surface || and transport properties of Luttinger liquids 
|J. In addition to extensive quantities which can be computed by means of the Bethe 
Ansatz |10| [III] , it is also interesting to analyse the behaviour of local observables in the 
presence of boundary conditions. In this paper we discuss a simple example of QFT with 
boundary conditions which has the advantage of being suitable for an exact analysis. The 
example we consider is the computation of one-point functions of the local operators of the 
Ising model defined in an half-space. Original lattice derivation of these quantities may be 
found in the refs. |12|] [|13 1 |l4j . Before proceeding in our computation, let us briefly recall 
some general aspects of the problem. 

The best way to compute correlation functions of a QFT with boundary is to take full 
advantage of the solution of the theory in the bulk and of the knowledge of the boundary 
state |§ |15[ . Namely, instead of considering the boundary placed at the spatial coordinate 
x = 0, it is most convenient to consider the boundary placed at t = and described by 
the boundary state | B > (Figure 1). 
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Figure 1. Boundary geometry for the computation of correlators. 
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Then, the correlation functions can be expressed as 

^ /n t + \ /n ( + w < \T t [Oi(xi,ti) ... ^ „ y ., n { ^ ^ 
< 1 {x 1 ,t 1 ) . . .O n (x n ,t n ) > = . (1.1) 
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In this geometry, the Hilbert space of the theory is as in the bulk and therefore, even in 
the presence of the boundary line, the local operators Oi can be completely characterized 
by the bulk Form Factors < (3\, . . ., f3 n \ Oi \ (3 n +i, ■ ■ ■ , Pm >, where fti are the rapidity 
variables, related to the two-dimensional momenta by p^ = rrii cosh/3^, p^ = mjsinh/3i. 

Let us focalize our analysis to the Ising model with a boundary in the low temper- 
ature phase (T < T c ) Q. The system can be described in terms of massive fermionic 
operators A((3) and A*(0) with the usual anti-commutation relations. The mass m is a 
linear measurement of the deviation of the temperature with respect to the critical one, 
m = 2n(T c — T). The spins on the boundary can be subjected to three possible boundary 
conditions, namely (a) they can be frozen to one of the fixed values ±1 ("fixed boundary 
condition"), (b) they can be completely free to fluctuate ("free boundary condition") or 
(c) they can be coupled to a boundary magnetic field h ("magnetic boundary condition" )0. 
The detailed analysis of these three possibilities is given in the original reference 0, and 
here we simply recall the results for the reflection S'-matrix R(f3) associated to each of 
them. For the free boundary condition we have R+(f3) = i coth — ^\ , for the fixed 

boundary condition R-((3) = — ztanh — and finally for the magnetic boundary 
condition 

^ , n \ ■ ■, ( & k — i sinh 3 ,„ „. 

Bh ( /3 ) = -,tanh(2-_j__ i J, (1.2) 

where k = 1 — Varying h, we can interpolate between the free boundary condition 
(h = 0) and the fixed one (h — > oo). 

Given an amplitude R(3), the corresponding boundary state can be then expressed 

asi 



B > = exp 



f-R{3)A\-3)A\3) 

Z7T 



> , (1.3) 



where R(3) = R — 0)- The structure of the boundary state \B) is such that the states 
entering its definition consist of pairs of particles with opposite rapidity (Cooper pairs). 



1 In order to simplify the notation, from now on we denote the free b.c. with the index (+), 
the fixed b.c. with the index ( — ) and the magnetic b.c. with (h). 

2 The most general expression of the boundary state may also involve an additional term 
relative to the zero mode. Here we discard it since it does not enter our following computations. 
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We can now use the Form Factors of the Ising model determined in [16] [I7| and the 
expression of the boundary state ( |1.3| ) to compute the correlation functions in the presence 
of the boundary. 



2. One-point and two-point functions of the energy operator 

The simplest correlation function is the one-point function of the energy operator 
e(x,t) which can be computed through the formula 

oo 

e (t) = J2<0\ e(x, t) | n >< n \ B > . (2.1) 

n=0 

The energy operator couples the vacuum only to the two particle state and for its matrix 
element in the euclidean space we have 

< | e(x, t) | Pi, p2 > = —2iimi sinh — — — 

x exp [—mt (cosh f3\ + cosh ^2) + ^mx (sinh j3\ + sinh (3 2 )} ■ 

(2.2) 

Hence the sum ( |2.1| ) consists of only one term and the one-point function of the energy 
operator can be expressed as 

POO 

e (t) = -im dPsmhpR(f3)e- 2mtcosh P . (2.3) 
Jo 

Its graphical representation is given in Figure 2. 

The one-point function does not depend on x, as a consequence of the translation 
invariance along this axis. The above integral reduces to closed expressions in terms of 
the modified Bessel functions in the case of free and fixed boundary conditions. Those are 
given respectively by 

e ^(t) = Tm [K x (2mt) ± K (2mt)} . (2.4) 

In the short distance limit mt — > 0, they have the scaling form €q (t) ~ -f^i, m agreement 
with the Conformal Field Theory prediction |jl~8 |. In the large distance limit, they instead 
decay exponentially with an extra power term, which is different in the two cases, i.e. 



B 



Figure 2. One-point function of the energy operator. 



Let us turn our attention to the one-point function of the energy operator in the 
presence of a boundary magnetic field. This is given by 



eo(t, h) = —m 



l'^ dP(cOsh(3 - 1) C ° S ^ + 1 e -2mtco S hp (2 6) 

J cosh p — K 



By varying h, we may switch between the free and the fixed boundary conditions, realizing 
therefore a Renormalization Group flow on the boundary. At first sight, however, it seems 
impossible to interpolate between the two functions eg (£), which present quite a different 
behaviour. (Notice, in particular, that the function e \t) diverges as +00 at t — > 
whereas the other correlation function 6q (t) goes to —00 near the origin.) To solve this 
apparent paradox, we have to correctly take into account the interplay between the two 
variables t and h. To this aim, it is convenient to write down a differential equation satisfied 
by eo(£, h). Observe that eq. (|2.6| ) can be initially expressed as 

e ( r ,h) = -m d{3 da (cosh/3 - l)(cosh/3 + «) e -™osh/3-a(cosh/3-K) ? ^ 
Jo Jo 

(r = 2mt) and the resulting integral can be expressed in terms of modified Bessel function 

dr]e Kri 



as 



n 
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Taking the derivative with respect to r in both sides of this equation, we obtain the 
following differential equation satisfied by eo(r, h) 



de (r,h) 

— — h Ke (r, h)=m 



Ki (r) 



[2.9) 



For h = 0, using the differential equation satisfied by the Bessel functions, we recover for 
eo(r, h = 0) the previous result eg (t), whereas e \t) is obtained by taking directly the 
limit k — > — oo in the differential equation. 

The differential equation is quite useful to analysing the short and large distance 
behaviour of the one-point function eo(r, h). For r — > 0, we can parameterize the solution 
as eo(r, h) ~ Cr~ s . Substituting into ( |2.9| ), it is easy to see that for any finite k, the term 
which determines the short distance behaviour in the RHS is m Kl , therefore C = —m 

r > 

and 8 = 1, independent of the value of the magnetic field. Hence, as far as h is finite, 
all the curves eo(r, h) follow at a short distance scale the behaviour associated to the free 
boundary condition. 

For r — > oo, we look for a solution of the form eo(r, h) ~ e~ r r~ {JJ J2^ =0 cik r ~ k ■ Substi- 



tuting into (|2.9| ) , using the large distance expansion of the Bessel functions, and comparing 
the power series we have then the following cases: 

(i) for k = 1, do = —2rriy^ and u = |, so that eo(r, h = 0) ~ aoe~ r r _1//2 . 

(ii) for — 1 < k < 1, we have ao = — m vlFx3/t anc ^ w = I' so ^^ a ^ e o /l ' > ( r ) ~ 
aoe -r r -3 / 2 . In this range ao < and the function approaches the real axis from be- 
low. 

(iii) for k = — 1 the first leading term of the expansion vanishes and for the next 
leading term we have a\ = — m^w'lj, u = | and the corresponding solution goes to zero 



16 V 2 ' ~ 2 

much faster, eo(r, /i = 2-y/m) ~ aie _r r -5 / 2 . 

(iv) for k < —I, we have ao = ~ m \f\^[^ j -> u = § an d ~ aoe _r r -3 / 2 . In this 
range ao > and the function approaches the real axis from above. 

The one-point functions relative to some of the above cases are shown in Figure 3. 

The result of this analysis is that for the boundary magnetic field smaller than the 
critical value h c = 2y / m, the one-point function does not have any zero at finite value of 
r whereas, for h > h c , it crosses the horizontal axes at a finite value of r, and its large r 
behaviour closely follows the behaviour relative to fixed boundary condition. 

The above picture is consistent with a Renormalization Group analysis. In fact, the 
boundary magnetic field is a relevant operator which therefore cannot affect the behaviour 
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Figure 3. One-point function of the energy operator versus (mt) for different 
values of the boundary magnetic field. Upper full line: fixed b.c. Lower full 
line: free b.c. Long dashed line: h < h c . Short dashed lines: curves with 
h > h c . 

of the correlation functions near the boundary, as far as it assumes finite values. Hence, 
sufficiently close to r = the boundary always appears subject to free boundary condition. 
The boundary magnetic field, however, influences the observables at large distance scales 
and when h > h c , the boundary always appears as subject to the fixed boundary condition 
for an observer placed at r — > oo. In an intermediate scale, there is a non-trivial crossover 
between the two different behaviours, which becomes rather sharp with increasing h. The 
function r (h), implicitly defined as the zeros of the one-point function eo(r(h),h) = 0, 
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Figure 4. Zeros of the one-point function €o(r(h), h). 



may be interpreted as a phase diagram of the theory which divides the free and the fixed 
boundary condition regions (Fig. 4). 

It is interesting to study the massless limit of the one-point function (|2.6| ). This 
corresponds to the physical situation of a critical Ising model in the bulk but in the 
presence of a boundary condition which breaks the conformal invariance of the model. 
With the change of variable y = 2mt cosh 0, the integral (|2.6|) becomes 



dye' 



y — 2mt y + 2mtK 



(2.10) 



'2mt V V + 2mt y - 2mtK 

In the limit m — > 0, we have 2mtK — > — h 2 t = —z and the one-point function of the energy 
operator can be expressed as 

1 



e Q (t, h) 



2t 



[l + 2ze z Ei(-z) 



(2.11) 



where Ei(z) is the exponential integral function. 

Let us consider now the connected two-point function of the energy operator in pres- 
ence of boundary magnetic field, 



G c (xi, t\\ x 2 , t 2 ; h) =< e(xi,ti)e(x 2 , t 2 ) > -e (ti, h) e (t 2 , h) 



[2.12) 
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By translation invariance along the x axis, this function depends on the difference x = 
X2 — x\. On the other hand, the presence of the boundary in the time direction implies 
a dependence of G on the two variables t\ and t 2 . It is convenient to define t = t 2 — t±, 



t == t 2 + 1± and r = \/ x 2 + t 2 . The correlation function ( |2.12| ) can be easily computed by 
using eqs. and the form factor (|2.2[) . Its closed form is given by 



G c (xi,ti;x2,t 2 ;h) = m 



2 



d „ , s „, _A 2 ( d 2 



K (mr) + F(mx,mt) + K (mr) 
o[mx) J \o(mt) 



d -\ f d - 
-=-F(mx,mt) — I Ko(mr) + — -F(mx,mt) 



d(mt) J \ d(mx) 



where we have introduced the function 

"+00 



(2.13) 



. % f , coshg + K q . . , nil . 

r{mp.mr) = — / aq — tanh - exp(— mr coshg + impsmhq) . (2.141 

2 J_ 00 coshg — k 2 

In the massless limit, the two-point function can be expressed as 

G c (xx, ti, x 2 , t 2 ; h) = 4j + t; h) F{—x, t; h) (2.15) 

where 

f°° 2q - h 2 

F(p,T',h)= dq- — — o exp(-rq + ipq) . (2.16) 
J 2q + h 2 

It is easy to see that for free and fixed boundary conditions, eq. (|2.15|) correctly reduces 
to the expression obtained by methods of boundary conformal field theory . 



3. One-point function of the magnetization operator 

Let us consider now the one-point correlation functions of the Ising magnetization 
field a(x, t). This can be computed by using the formula 

00 

<r (t) = J2(0\a{x,t)\n)(n\B) , (3.1) 

n=0 

and by translation invariance, it only depends on t. We denote by <jQ + \t) the one-point 
function relative to the free b.c, by Op (t) the one-point function with fixed b.c, and 
finally by oo(t, h) the one-point function in the presence of an arbitrary boundary magnetic 
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field h. In order to compute these one-point functions, we need the bulk Form Factors of 
the magnetization operator I, given by [ I(| Jl7] : 



<0|a(0)|ft ■■■&„>= z n n tanh (^y^) • (3 - 2) 

Since the boundary state \B) consists of pairs of particles with opposite rapidity, we have 
to specialize the previous formula and evaluate 

) = JJtanhA J] ^tanh^-^tanh^-^V = 

i=l i<7 ^ ' 



(0|o-(0)| -PuPx, 



nr n tanh ^ n 



/ cosh /3i — cosh /5j x 



V cosh A + cosh Ba 

t=i t<j v ri ^ J 

(3.3) 

Noting that 



n( cosh A - cosh Ba \ 2 , /2a/ cosh cosh \ TTr//n ^. . . 

— r-T r^- = det v , n — — ^ = det W(3i,pA , 3.4) 
. V cosh A + cosh A / \ cosh A + cosh Pj 1 

-(±), 



we can express the one-point function Uq (£) as a Fredholm determinant, namely, 

-+oo jo /-+00 



o ±} W=EA/ ^ IT'" / "f^fn^tanh^^^))) detW(A,A 
= Det(l-*±V( ±) (t)) , 



(3.5) 

where the kernels are given by 

= 'SI'S) • e<±,(ft<) = ^M=Te— -0. (3.6) 

In the above formula, the physical correlator is obtained for z± = ±1/(2tt), but it is 
convenient to regard the above expressions as functions of the parameters z±. 

The one-point function 0"o(£, h) can also be cast into a Fredholm determinant form as 

a (t,h) = Det(l -zV(t,h)) , (3.7) 



3 With our normalization of the Form Factors, the conformal limit of the two-point func- 
tion in the bulk is given by (a(r)a(O)) = F 2 r~ 1/A , where T = 2~ 1/12 e 1/8 A~ 3/2 m~ 1/8 = 
0.73642778 m~ 1/8 , with A = 1.282427 the Glasher constant. 
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where the kernel is now given by 



V(Pi,f3j,t,h) 
E((3,t,h) 



E(Pi,t,h)E(Pj,t,h) 
cosh Pi + cosh Pj 

( cosh _ 1) ^shP + K e _ mt cosh p 

cosh P — K 



(3. 



In this case, the physical correlator is obtained for z — l/(2ix). 

It is quite easy to show that the crossover phenomenon previously discussed for the 
one-point function of the energy operator also takes place for o" (i, h). Let us initially 



consider the short distance behaviour of eqs. (|3.5| ) and ( |3.7|) . As proved in appendix A, in 
the limit mt — ► 0, we have 

Concerning the behaviour of the one-point function ao(t, h) in the limit mt —* 0, its power 
law singularity is independent of h, as far as h is finite, and coincides with that one relative 
to the free boundary condition, i.e. ao(t,h) ~ (2t) 3 / 8 . Hence, varying h, all the curves 
er (t, h) follow at short distance the power-law behaviour dictated by the free boundary 
conditions. 

Let us evaluate now the large distance behaviour. For the fixed boundary conditions, 
we have 



i + 



i 



-2mt 



+ 



(3.10) 



16v^F (mt) 3 / 2 

i.e. this function approaches the bulk expectation value from above. For the free boundary 
condition, we have instead 



1 



1 



-2mt 



2 v ^F(mt) 1 /2 



+ 



(3.11) 



i.e. the bulk expectation value is reached from below. In the presence of a boundary 
magnetic field, the asymptotical approach to the bulk vacuum expectation value is given 
for ft ^ — 1 by the expression 



whereas for k 



a Q (t, h) 



a (t, h) 



1 



1 1 + k e 



-2mt 



16^F 1 - k (mt) 3 / 2 



+ 



-2mt 



1 - 



128v^F (mt) 5 / 2 
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+ 



(3.12) 



(3.13) 







0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 



Figure 5. One-point function of the magnetization operator versus (rnt). Upper 
full line: fixed b.c.. Lower full line: free b.c Dashed line: curve with h > h c . 

From ( |3.12|) we see that for — 1 < k < 1, the curve reaches its asymptotic behaviour 
from below, whereas for k < — 1 from above. At k = — 1 there is a faster decreasing 
behaviour. Hence, also for the one-point function of the magnetization operator there is 
the typical Renormalization Group crossover between free and fixed boundary conditions, 
moving away from the boundary (Fig. 5). 

4. Differential equation for the magnetization one-point function 

Bariev has shown that for the free and fixed boundary condition, the one-point func- 
tions of the magnetization can be expressed in terms of a solution to the Painieve difler- 
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ential equation [13 . He started from a lattice description, from which he obtained the 
same Fredholm determinant expressions as we directly obtained by using the form factor 
approach. In this section we derive differential equations for arbitrary magnetic field h 



using the general techniques developed in |n| |2(J • As we will show, these will turn out to 
be coupled non-linear partial differential equations in the variables t and k. 

Making the change of variable u = cosh /3, we can initially write the one-point function 
<Jo(t, h) as 

°° 1 f + °° din f +oa dn 
(0\a(t)\B(h)) = det(l + V) = J2- / " / %T detV(«i,^,t) (4.1) 

n=0 ^ 1 ** X 



where 



V{ui,Uj,t) 



e(u, t) 



e(ui,t) e(uj,t) 



Ui + Uj 



u — 1 



1/4 



u + 1 J v /<' — a 
One first needs the resolvent R±, defined as follows: 

(l-i2±)(l±V) = l • 



K + u e -mtu 



(4.2) 



(4.3) 



In the following, it will always be implied that Fredholm operators are multiplied in the 
usual way, e.g. 



{RV)(u,v) 



+00 



dw 
2^ 



R(u, w)V(w, v) 



and similarly for multiplication of operators with functions: 

"+00 



(Re)(u) = J —R(u,v)e(v) . 



(4.4) 



(4.5) 



Resolvents for the kernels of the type ( |3.6| ) were described in [21]. Define functions / ± (w) 
by the formulas 

(1 ± V)/± = e ^ /± = (1 - R ± )e . (4.6) 



Then 



where 



H(u,v) 
F{u,v) 



R±(u, v) = H(u, v) ± F(u, v) , 

f-(u)f+(v)-f+(u)f-(v) 



2(u-v) 
f + (u)f-(v)+f-(u)f+(v) 
2(u + v) 
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(4.7) 



(4- 



The next step is to identify variables with the property that derivatives of the kernel 
V with respect to such variables acts as a projector, i.e. removes the l/(u + v) factor. The 
spacial variable t obviously has this property: 



dtV(u, v) = — me(u)e(v) . (4.9) 

One also has 

d K V(u, v) = (uv - K 2 )e'{u)e'{v) , (4.10) 

where 

e'(u) = ■ (4-H) 

For future convenience, define also 

e(u) = -^- 2 e{u) . (4.12) 
The rest of the derivation of correlation functions parallels closely the treatment of 



Ising 2-point functions given in |21J] with z, z replaced by r, k. Let 

r± = Det(l± V) . (4.13) 

Then 

d t \ogr± = ^Trlog(l ± V) = Tr(l - R±)d t V 

(4.14) 

= =Fm P± , 



where 



P ± =<e,/ ± >^ / ^eH/^w) . (4.15) 



Next, we compute d K P±, for which we need d K f^. Taking the derivative of ( |4.6| ), one 
obtains 

(l±V)dJ ± ±d K Vf ± = e . (4.16) 

The latter implies 

(1 ± V)dJ ± =F K 2 e' < e', f ± >±e< e, /± >= e . (4.17) 

Define the additional functions 

/± = (1 - R±)e , f f± = (1 - R±)e' . (4.18) 
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Then multiplying ( fll7|) by (1 — R±), one obtains 

dJ ± = f ± ±K 2 D ± f' ± T D ± f ± 



where 



One then finds 



D ± =< e'J± >=< e,/ ,± > 
D ± =<eJ ± >=< e,f ± > 



(4.19) 



(4.20) 



(4.21) 



d K p ± =<d K ej ± > + <e,dj ± > 

= 2D ± ±k 2 (D ± ) 2 t(D ± ) 2 . 
We can also obtain expressions for d^D^ and dtD^ 1 . For this, one needs dtf^- This 



is the same as in [^T 



d t f ± =m [±(P+ + p-)f ± -ufT] 



One then has 



d t D ± =<d t e',f ± > + <e',d t f ± > 



= m 



±(P + + P~)D ± — D + — D 



± _ n+ 



(4.22) 



(4.23) 



Similarly, 



d f D ± = m 



±(P + + P~)D ± - (P+ + P~) - k 2 (D+ + D~) 



(4.24) 



Finally, there is one non-trivial constraint among the D, D potentials, which is proven 
in the appendix B: 

D+ = D~ - D+D~ + k 2 D+D~ . 



(4.25) 



We can summarize the above results as follows. Define 



Then, 



P = p+ + p~ , 

C = 2 + £>- -D + , 

D = D + + D~ , 

d t log (t + /t_) = -mP 

d t log (t + t_) = -mQ , 

d K P = k 2 DD - CC , 
1 



Q = P+ - P~ , 
C = -(D + + D- 
D = D + -D~ . 



(4.26) 



(4.27) 



d K Q = 2--\C* + C A ) + —[D z + D< 



K 
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and 



d t D = m 2C + PD 
d t D = mPD , 
d t C = mPC , 
d t C = m [PC + 2k 2 D] 



(4.28) 



The constraint reads 



C 2 -C 2 + n 2 {D' 



4 . 



(4.29) 



The equations ( |4.27| ) indicate that dtd K log(r±) can be expressed in terms of the four 
functions C, C, D, D. In turn, the function P can be eliminated from Q4.28|) to obtain 
differential equations involving C,C, D, D only. For instance, expressing P = dtD/mD, 
and using the expression for d K P, one finds: 



Dd K d t D - d t Dd n D = mD (k DD - CC). 



(4.30) 



Three additional differential equations of the same kind can be similarly obtained by 
eliminating P from each of the equations in (4.28). 



5. Concluding Remarks 

In this paper we have studied the effects induced by a boundary Renormalization 
Group flow on the correlation functions of the two-dimensional Ising model. Using the 
form factors of the local energy and magnetization operators, we have computed their one- 
point functions in the presence of a boundary magnetic field as sums on the intermediate 
states of the bulk theory. The integral representations of the correlation functions can be 
used to derive the differential equations which they satisfy. In view of the "free" nature 
of the energy operator in the Ising model, its one-point functions in the presence of a 
boundary magnetic field is a solution of a linear non-homogeneous differential equation, 
given in ( |2.9| ). On the contrary, the magnetization operator is an interacting field and for 
its one-point function in the presence of a boundary magnetic field, we obtain the system of 
non-linear differential equations (|4.27|) , (f4.28| ) ,with the constraint expressed by (|4.29| ) . The 
rather different mathematical structure of these differential equations reflects the distinct 
dynamics of those fields. 

Some simplications in the equations are expected in the massless limit of the theory. 
This corresponds to the Ising model with an infinite correlation length in the bulk but 
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with a violation of the conformal symmetry induced by the boundary interaction. We 
have analysed this interesting situation for the cases of the one-point and the two-point 
functions of the energy operator. Concerning the one-point function of the magnetization 
operator, it is known that it satisfies in the massless limit a linear differential equation of 



by analysing the massless limit of the system of non-linear differential equations derived 
in this paper. 
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Appendix A. 

In this appendix we calculate the power law singularity of the one-point function of the 
magnetization operator. Let us start by analysing the free and fixed boundary conditions. 
An equivalent expression for ^^(t) is given by 



hypergeometric type It would be interesting to obtain directly the result of ref. [E3 



Acknowledgements 



a^\t) = det(l - z±V (±) ) = exp [Trlog(l - z±V (±) )] 



00 yti p+oo r+oo n 

n=l U J ° J ° i=l 



— 2mt cosh Pi 



cosh Pi + cosh P i+ i 



cosh Pi ± 1 



71=1 



(A.l) 



(where P n+1 



Pi). Using now 




da e 




(A.2) 



the coefficients Z± (t) in (|A.1|) can be expressed as 




(A.3) 



x [K 1 (a 2 + a 3 ) ± K (a 2 + a 3 )] • • • [K x (a n + ax) ± K (a n + an)] , 
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where Ki(x) and K (x) are modified Bessel functions. For mt — > 0, the leading singularity 
of these expressions comes from 



Z±\t)~ \ dai--- da n Ki(ai + a 2 )K 1 (a2 + a 3 ) ■ ■ ■ K 1 (a n + ai) . (A.4) 

J mt J mt 

It is easy to estimate the sum of the most singular terms. To do this, let us express a^ft) 
in terms of the eigenvalues X±(t) of the integral operators V± and their multiplicity ai(t) 
as 

^ ± )(t) = n(i-^>(t)f <l) . (A.5) 

i=l 

As far as (mt) is finite, the kernel is square integrable. However, when (mt) — > 0, the 
operator becomes unbounded. In this limit, the eigenvalues becomes dense in the interval 
(0, oo) according to the distribution 

Hp) = , (A.6) 

cosh 7rp 



iln^. 

71" mx 



whereas, from Mercer's theorem, their multiplicity grows logarithmically as ~ — In 
Therefore 

loga ( ±} (t) =x±logt (A.7) 

where the critical exponents x± of the magnetization operator relative to free and fixed 
boundary conditions are given by 

1 f°° , , ( \ 2nz ± \ 1 1 2 



x± = / dp In 1 — = - - — rarccos (-2ttz±) . (A. 8) 

n z J \ coshp / 8 2n z 

Substituting the values of z± we obtain x = 3/8 for the free b.c. and x = —1/8 for the 
fixed b.c. 

Repeating the same analysis for the one-point function in the presence of the boundary 
magnetic field h by using the formula (|A.1|) , it is easy to see that, as far as h is finite, the 
most leading singularity of the corresponding is always given by ( |A.4j ) and can be 
estimated according to the equation ( |A.8| ) with z = 1/(2tt), i.e. x = 3/8. 
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Appendix B. 

Here, we relate the functions f^, / /=t: to the functions / ± and prove the constraint 

(S). 

Using the identity 



11 1/1 1 



u + v v + k u — k \v + k u + v 
one first easily shows that 



(B.l) 



From the definition, 



V 2 — K 2 



Factorizing v 2 — k 2 = (v + k)(v — k) in the above equation, and applying ( |B.2| ) twice, one 
obtains 

= - J-j [(1 ± Q*)/* =F «D T / ± =F uf^D*] , (B.4) 



where 



Using the identities 



± =< ± _e(u)_ >= 5± + KjD± (B 5) 

"U — K 



H{u, v)v = uH{u, v)+ l - (/+(«)/"(«) - /-(«)/+(«)) , 
F(u,v)u = -«F(«,v) + ^ (f+(u)f-(v) + f-(u)f+(v)) 



(B.6) 



one also finds 

f±(u)=uf'*TD*f ± . (B.7) 

The constraint (|4.25|) follows from using ( B.7|) , ( B.4j ) to express J 1 * 1 in terms of J 1 * 1 in 
the equation D ± =< e, f ± >=< e, f ± >. 
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